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Abstract 

High resolution nanoscale imaging in liquid environments is crucial for studying molecular 
interactions in biological and chemical systems. In particular, electron microscopy is the gold- 
standard tool for nanoscale imaging, but its high-vacuum requirements make application to in¬ 
liquid samples extremely challenging. Here we present a new graphene based wet cell device 
where high resolution SEM (scanning electron microscope) and Energy Dispersive X-rays 
(EDX) analysis can be performed directly inside a liquid environment. Graphene is an ideal 
membrane material as its high transparancy, conductivity and mechanical strength can support 
the high vacuum and grounding requirements of a SEM while enabling maximal resolution and 
signal. In particular, we obtain high resolution (< 5 nm) SEM video images of nanoparticles 
undergoing brownian motion inside the graphene wet cell and EDX analysis of nanoparticle 
composition in the liquid enviomment. Our obtained resolution surpasses current conventional 
silicon nitride devices imaged in both SEM and TEM under much higher electron doses. 
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Introduction 


Imaging in a liquid environment is important across a wide range of research fields from physics 
to biology. The nanoscale imaging of such systems drives new insights in molecular and biologi¬ 
cal theory.^ New experiments enabled by wet-cell technology include live imaging of antibodies 
and bacteria to understand immune response, and in-situ imaging of crystals to understand growth 
kinetics.13 Most of the imaging is performed with electron microscopy such as SEM (Scanning 
Electron Microscope) or TEM (Transmission Electron Microscope)In particular, SEMs are 
widely available and accessible to most researchers. While SEMs offer quick and high resolu¬ 
tion nanoscale (2-lOnm) imaging, the high vacuum operation conditions (<10 ^ Torr) of these 
instruments make the imaging of liquid environments challenging .1^ Systems that operate at high 
pressures such as environmental SEMs (E-SEMs) are specialised tools requiring the use of water 
vapour to purge and replace air in the specimen chamber. Moreover, the electron beam in such 
systems scatters from the introduced vapour resulting in limited resolution.!^ 

Conventional wet cells are based on sealing liquid samples behind a 30-150 nm silicon nitride 
window.!^^^ While this approach has proved effective, the resolution is fundamentally limited by 
the necessity of using relatively thick nitride membranes. Experiments have obtained a resolution 
of only around 20 nm for a membrane thickness of 50 nm in an SEM.^^The fabrication of thinner 
nitride windows with thickness below 50 nm is challenging, requiring special techniques to con¬ 
trol the etching rate and achieve etching uniformity.!^ As the nitride membrane becomes thinner, 
the windows become too fragile to handle. Silicon nitride wafers are also electrically insulating, 
requiring the sputtering of a thin layer of conductive material such as gold for electrical leads or to 
ground the sample.!^ Ultimately, nitride based windows cannot be extended to thicknesses below 
a few nanometres. This is a very crucial technical limitation, limiting not just resolution but signal. 
Eor example, the need for relatively thick nitride windows obviates application of standard SEM 
techniques such as Energy Dispersive X-Ray (EDX) due to the absorption of signal by the thick 
membrane. 

Here we present a graphene wet cell for SEM imaging under a high vacuum environment. 
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Graphene is an atomieally thiek layer of earbon atoms (0.34 nm thiokness)l^with exeeptional prop¬ 
erties ineluding high meehanieal strength, high thermal and high eleetrieal eonduetivity. Graphene’s 
atomie thiekness makes the material an optimal imaging window enabling maximum resolution 
and signal. In partieular, graphene allows for the eolleetion of low energy seeondary eleetrons 
as opposed to just baekseattered eleetrons performed in most SiN wet eell imaging studies.^ 
This greatly improves the signal and resolution of the images.aAI Graphene’s meehanieal strength 
prevents breakage of ~5 jum membranes under vaeeum eonditions. Graphene’s high thermal eon¬ 
duetivity allows exeess heat generated from the beam to dissipate quiekly without damaging the 
sample. Finally, graphene’s high eleetrieal eonduetivity obviates the need for an additional metal 
eoating for grounding. The graphene membrane also provides eonvenient eleetrieal leads for 
voltage and eurrent inputs for adding eleetrieal bias in experiments. Previous groups have used 
graphene oxide membranes for imaging.^ However it is ehallenging to eontrol the homogeneity 
in the graphene oxide membrane aeross the window and, at around 20 nm thiek, they are eompa- 
rable in thiekness to nitride. Using ehemieal vapour deposition (CVD) with earefully eontrolled 
growth eonditions we ean ensure that there is a single layer graphene membrane 

Our single-layer graphene wet eell deviee enables dynamie imaging in a SEM. In partieular, 
we observe Brownian dynamies of Au-NP’s transiently binding and unbinding at the surfaee of the 
graphene. While Brownian motion of Au-NP’s has been observed previously in a TEM using a 
graphene sandwieh assay, developing a moleeular in-liqiuid imaging eapability in an SEM has key 
praetieal and fundamental benefits .^^SEM’s are more available, eheaper and more versatile tools 
that permit introduetion of mueh larger samples. Eor example, large (1-10 em size) miero/nano flu- 
idie deviees eould be easily introdueed into an SEM and wet eell imaging eould then be performed 
as part of routine deviee operation. In partieular, as there is no eonstraint on sample thiekness in an 
SEM, an SEM-based wet eell ean ineorporate mueh deeper fluidie ehannels without loss of signal, 
signifieantly simplifying wet-eell design. Moreover, additional sample material ean be potentially 
pulled in from deeper in the eell. Eor example, we show that eontinuous seanning attraets Au-NP’s 
to the graphene interfaee. Einally, SEMs are outfitted with a wide range of surfaee eharaeterisation 
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tools (for example, EDX). We show that, using our graphene wet eell deviee, these tools ean then 
be adapted to study the wet eell environment. As an example, we able to obtain an EDX speetrum 
of Au NPs in liquid. 


Sample preparation 


Our fabrieation proeess is divided into three steps, the fabrieation of the silieon nitride substrate, 
the growth and transfer of the graphene and the wetting and sealing of the deviee for SEM imaging. 
An illustration of the deviee is shown in Eigure 


a) 



Silicon substrate 


• Gold nanoparticles 
Liquid Reservoir 


Eigure 1: a) A schematic of our graphene wet cell device, b) SEM image of the liquid environment imaged 
through the graphene membrane micropore. The graphene is positioned on top of a circular aperture etched 
through the SiN membrane, c) Schematic of device as viewed from the side. The liquid sample, held in the 
400 /rm fluid reservoir sandwiched between the graphene membrane and kapton tape, consists of deionized 
water with Au nano particles. The figure is not drawn to scale. 
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The first step is substrate fabrieation. Our substrate is a 400/im thiek (110) silieon wafer eoated 
with a 180nm thiek nitride membrane and divided into 2x2mm dies. The wafer was patterned 
with photolithography and etehed in KOH from the baek to produee a 70x70 jum residual nitride 
membrane in the middle of eaeh die. The KOH etehed apertaure also serves as a reservoir for the 
liquid sample. Lastly, a 2 jum diameter hole was etehed through the middle of the free standing 
nitride membrane to form the graphene viewing window. 

Graphene was grown using Chemieal Vapour Deposition (CVD) on a 25 jum thiek eopper foil 
with a growth temperature of 1050 ° C at a pressure of lOOmTorr and a flow of 4 seem of CH 4 .^ 
Our eustom-built CVD system is based on a vertieal furnaee. Two gas tubes feeds into the top of 
a 2.5 em wide vertieal quartz tube to provide the flow of gases. The quartz tube is lowered into 
the oven during the growth and the growth time is approximately 1 hour. Before the growth, the 
eopper foil was first annealed in a flow of 12 seem of hydrogen for an hour to strip the oxide layer 
on the foil. The CVD synthesized graphene was then spin eoated with a thin supporting layer 
of polymethylmethaerylate (PMMA) layer and the Cu substrate was etehed away in a solution 
of 0.1 M ammonium persulfate ((NH4)2S208). The sample was transferred by inserting a glass 
slide into the ammonium persulfate solution, using the slide to seoop out out the freely floating 
graphene membrane and depositing the graphene bearing slide into a beaker of de-ionised water. 
To eompletely remove the ammonium persulfate residues, the sample was transferred into another 
elean beaker of de-ionised water before being transferred onto the top side of the silieon nitride 
wafer sample to eover the 2 jum holes. Graphene produeed using the same growth eonditions was 
transferred onto Si02 wafers for Raman speetroseopy to eonfirm that the graphene was indeed 
monolayer. 

Finally, the sample was ready to be wetted and sealed. Gold nanopartieles (Au-NP’s) 20 and 
50 nm in diameter were used to eharaeterise the fluid eell. We ehose Au partieles as they are eom- 
mereially available in a wide variety of sizes and ean potentially be used as eonduetive biologieal 
labels.^ The Au-NP’s were diluted 1:20 from stoek solution in DI and then the nanopartiele eon- 
taining solution was degassed for an hour. Degassing was erueial to ensure proper wetting and 
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to decrease the formation of gas bubbles during imaging. After degassing, several microliters of 
solution was pipetted into the reservoirs and the wafer sample was sealed with Kapton tape on the 
back side. The device was then rinsed in acetone and isopropanol to dissolve the PMMA support¬ 
ing layer on the graphene. The imaging of the device was then done using a FEI-F-50 SEM in 
the standard high vacuum mode at 10 ^ Torr using a secondary electron detector of the Everhart- 
Thomley type. The graphene membrane remained intact at this operating pressure of 2.2 x 10^^ 
Torr. The primary electron energy used for imaging is 10 KeV. Under these imaging conditions, 
the escape depth of secondary electrons in water should be in the order of 10s of nm.l^ 


Observation of Nanoparticle Dynamics 

While many Au-NPaAZs are non-specifically bound to the membrane, we observe Brownian dy¬ 
namics of Au-NPaAZs floating in solution below and undergoing transient interactions with the 
membrane. These dynamics are recorded over several minutes using a screen capture program. 
Eigurej^gives an example of bead motion. Beads are observed to be diffusing in and out of contact 
with the membrane surface, confirming that they are indeed contained in a liquid environment. The 
particle trajectories are recorded using a custom tracking program.^ 

In the absence of confinement, the gold nanoparticles are expected to undergo Brownian motion 
in water, characterized by a diffusion constant: 


D 


hT 

6lir\r 


( 1 ) 


where k^, is the Boltzmann’s constant, T is the temperature (300 K), r\ is the viscosity of DI 
water (lx 10^^ Pa- S), r is the radius of the beads (25 nm). Eor our image frame rate of l/t5=29Hz, 
this leads to a corresponding mean diffusion length of Ljj = y/Dfs — 550 nm at room temperature 
(300 K). Hence, within one image frame the particles are expected to approximately hop 1/4 of the 
length of the nanopore. Eigure suggests that we indeed see fluctuations on that scale. However 
we also observe two additional types of behaviours. Particles can be permanently bound to the 
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Figure 2: Image time-series showing Au-NP dynamics in our graphene wet-cell device. Beads are observed 
to be diffusing in and out of contact of the graphene nanopore. The white bar indicates 500 nm. 

membrane over the course of the imaging time and can also diffuse in and out of contact with the 
membrane, interacting transiently with what appears to be “sticky sites". This sticking behaviour 
can be quantified by a plot of occupation probability p{x,y). The occupation probability is taken 
by integrating the total number of frames a bead appears at a certain location normalised over the 
total number of frames of the video. Figurej^shows the occupation probability for the same device 
with a spatial resolution of 10 nm and time resolution of 25 ms, clearly indicating the existence 
of strong trapping sites that permanently bind beads and weaker trapping sites that give rise to 
transient interactions. 
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Figure 3: a) Plot of the integrated (over 15s) normalized nanoparticle occupation probability across 
the graphene membrane. The nanoparticle occupation is defined as a brightness of 70 % or more. 
Some spots show an occupation of unity, meaning that beads are bound to the membrane at these 
positions for the entire duration of the movie, b) Zoomed image of the upper left comer with arrows 
indicating the positions of the time traces in (c). Each pixel shown corresponds to an integrated 
area of 30x30nm^ at a frame rate of 29Hz. The scale bar denotes 500 nm in length. 


The non-uniformity of the occupation probability suggests that the graphene membrane varies 
with regards to its physical and chemical reactivity towards nanoparticles. One possible source of 
non-unformity are the existence of grain-boundaries in the graphene layer .1^ These grain bound¬ 
aries are imperfections in the graphene lattice due to differently orientated growth directions. The 
grain boundaries from previous studies are spaced roughly the same distance apart (~ 1/im) as 
the observed sticky sites. Another possible source of non-uniformity is the presence of graphene 
“wrinkles" arising from the growth conditions on the inhomogeneous surface of the copper foils.^ 
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The wrinkles form valleys in the graphene sheet allowing beads to be drawn in through attraction 
by van der Waals forces (which has also been observed in other wet cell applications).To reduce 
this effect, we repeated the experiment with PEG (polyethylene glycol) coated Au beads as shown 
in Figure Indeed we observed a suppression of the adhesion of Au particles to the graphene 
membrane with a reduction of the density of stuck beads upon imaging ^ 



Figure 4: a) The graphene pore before imaging in SE mode. The red boxes indicate the areas selected to 
be scanned for 2 minutes, b) The image of the pore after scanning, c) Image of a graphene membrane on a 
different device that has been repeatedly scanned by the electron beam. Note that the beads are drawn to the 
surface. The white scale bar indicates 1 /rm. 


Our dynamic SEM imaging capability, performed over a deep sample reservoir, enables us 
to demonstrate that continuous scanning draws beads to a scanned region of the membrane. We 
selected two areas (0.7x0.2/rm) on the graphene membrane, marked by the dashed red boxes in 
Figure]^ a). The area was then scanned continuously at 5 KeV for 2 minutes. Beads were observed 
to diffuse onto the graphene membrane in those areas (Figure |^b). Note that Figure |^b) is exactly 
at the same spot as Figure]^ a). The image looks different because the graphene membrane was 
observed to be deforming from the continuous scans. In Figure |^c), we scanned the entire pore on 
a different device for several minutes to draw beads onto the graphene membrane. This dynamic 
beam-induced attraction of the beads might arise from electrostatic charging of the membrane, 
possibly related to chemical modification of local impurities such as PMMA residues (resulting in 
charge trapping).!^ In addition, the PMMA from the supporting layer in the transfer process may 
not be totally removed during acetone-based dissolution process. Finally, space charge transiently 
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deposited in a nanoseale region beneath the graphene by either eleetron depositon or seeondary 
electron generation might induce polarization forces on the beads The stability of the graphene 
membranes greatly varies due to these effects. We observed cells that were stable for 2-15 minutes 
under 10 KeV. Future cells can be improved by optimizing the imaging conditions and reducing 
contamination of the graphene membrane. 


Energy Dispersive X-ray Spectroscopy 

The use of EDX in a wet cell could potentially allow for positive chemical identification of ele¬ 
ments in a liquid environment. EDX, however, cannot be performed in a standard silicon nitride 
wet cells due to the thickness of the nitride layer that absorbs emitted radiation. Here we show that 
graphene membranes enable EDX-based analysis in liquid enviornments. 

We attracted diffusing Au NP’s with the electron beam to the surface of the graphene membrane 
and performed EDX. Eigure shows spectrums taken at two different locations. The first EDX 
location (Eigure a)) is for a graphene bead located underneath the membrane layer. The second 
location (Eigurej^b)) is for a Au-NP resting on top of the graphene-silicon nitride wafer away from 
the membrane. To ensure that the EDX spectra correspond to beads in liquid and not on the surface 
of the membrane, we performed EDX only on beads that had freely diffused onto the surface of 
the graphene membrane during imaging and were stuck there during the EDX. Remarkably, we 
only see a 30% reduction in the integrated intensity for the Au signal under the graphene vis-a- 
vis the control spectrum. The source of the attenuation may be due to absorption of the signal by 
surrounding water or contamination deposited by the electron beam during the EDX measurement. 
Despite the attenuation, we are still able to positively identify the in-liquid particle composition. 
We also observe a much lower but non-zero silicon peak coming from location (a). The peak arises 
from the silicon background scattering from the hole edges. In addition, we observe a weak copper 
peak on the suspended graphene membrane, likely arising from copper used in the growth process 
that is not completely removed. We were able to perform multiple EDX measurements without 
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any degradation of the graphene membrane. 
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Figure 5: EDX speetra for a bead underneath the membrane (blue line) and at a eontrol loeation 
away form the membrane and on top of the nitride film (red line). Both speetra were integrated 
over 30 s and the bead-in-liquid speetrum is offset by 300 eounts for ease of eomparison. Inset 
a) Seeondary eleetron image of loeation of bead-in-liquid. Inset b) Seeondary eleetron image of 
eontrol bead. Both seale bars indieate 500 nm. 


Resolution 

To determine the resolution of the Au-NP’s under the graphene wet eell we imaged 20 nm beads 
bound to the membrane. The intensity line profile of eaeh Au-NP (see Figure was extraeted and 
the resolution was determined from the edge-width over whieh the Au-NP’s intensity rose from 
20% to 80% of its maximum height.l^ Averaging over five beads, we find the resolution of the 
Au-NP’s in our wet-eell to be 5 ± 3 nm (error is standard deviation on mean over beads measured). 
The eontrast to noise ratio (CNR) is defined as : 

CNR = — (2) 

where S is the peak signal and is the standard deviation of the baekground noise.^^ We 
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obtained a value of 7 ± 1. These results eonfirm that graphene leads to improved resolution: our 
resolution is higher than the 20 nm reported in the 50 nm silieon nitride membrane and eomparable 
to the resolution (~ 5 nm) obtained under mueh higher eleetron imaging eonditions sueh as at 200 
Kev in a TEM.I^ The high eontrast to noise (CNR) ratio also makes it possible for us to observe 
and reeord movements in the liquid environment. 



Pixel 

Figure 6: a) A secondary electron image of 20 nm Au particle non-specifically absorbed to the graphene 
window, b) A close-up image of a bead selected as an example, c) Intensity line profile. The edge width 
is determined from the 20-80% rise in intensity at the profile edge. A resolution of 5 ± 3 nm with a CNR 
of 7 ± 1 is obtained from averaging results over five beads. The white and black scale bars in a) and b) 
corespond to 500 nm and 25 nm respectively. 


Conclusion 

In conclusion, we have demonstrated that single layer CVD grown graphene is very promising for 
SEM based wet eell imaging, enabling dynamie imaging of Au-NP undergoing brownian motion 
in aqueous solution and EDX measurements in liquid. In partieular, our wet eell ean be used in a 
conventional SEM without the need for instrument modifieation. In the future, opposed to previous 
graphene sandwitch studies, our wet-eell ean be in principle adapted for nanofluidic experiments 
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with nanochannels etched in place of the fluid reservoirs. Such systems might enable nanoconfine¬ 
ment based single molecule manipulation combined with SEM imaging, giving rise to new types 
of single-molecule analytical devices based on electronic as opposed to optical imaging. Our EDX 
results are particularly significant in this context: one can envision future experiments that use 
biomarkers with differential chemical composition to tag a range of DNA modifications enabling 
more efficient multiplexing. 

Acknowledgement 

The authors acknowledge financial support from Eonds de recherche Nature et technologies (EQRNT). 
The authors thank the staff at the Eacility for Electron Microscopy Research (EEMR) at McGill 
University for the use of their electron microscopy instrument and technical support. 

References 

(1) Ackerley, C.; Nielsen, C.; Hawkins, C. Microscopy and Microanalysis 2006 , 12, 428-429. 

(2) Ishii, M.; Nakamura, H.; Nakano, H.; Tsukigase, A.; Harada, M. Langmuir 2005 , 21, 5367- 
5371. 

(3) Niu, W.; Zhang, E.; Xu, G. Nanoscale 2013 , 5, 3172-3181. 

(4) Kolmakov, A.; Dikin, D. A.; Cote, E. J.; Huang, J.; Abyaneh, M. K.; Amati, M.; Grego- 
ratti, E.; Gunther, S.; Kiskinova, M. Nature nanotechnology 2011 , 6, 651-657. 

(5) Tsuda, T.; Nemoto, N.; Kawakami, K.; Mochizuki, E.; Kishida, S.; Tajiri, T; Kushibiki, T.; 
Kuwabata, S. ChemBioChem 2011 , 12, 2547-2550. 

(6) Bogner, A.; Jouneau, R-H.; Thollet, G.; Basset, D.; Gauthier, C. Micron 2007 , 38, 390-401. 

(7) Muscariello, E.; Rosso, E; Marino, G.; Giordano, A.; Barbarisi, M.; Cafiero, G.; Barbarisi, A. 
Journal of cellular physiology 2005 , 205, 328-334. 


13 



(8) Peckys, D.; Werner, U.; de Jonge, N. Microscopy and Microanalysis 2013, 19, 138-139. 

(9) Donald, A. M. Nature materials 2003 , 2, 511-516. 

(10) Yang, L.; Zhu, Z.; Yu, X.-Y; Rodek, E.; Saraf, L.; Thevuthasan, T.; Cowin, J. P. Surface and 
Interface Analysis 2014 , 46, 224-228. 

(11) Mueller, C.; Harb, M.; Dwyer, J. R.; Miller, R. D. The Journal of Physical Chemistry Letters 
2013 , 4, 2339-2347. 

(12) Thiberge, S.; Nechushtan, A.; Sprinzak, D.; Gileadi, O.; Behar, V.; Zik, O.; Chowers, Y; 
Michaeli, S.; Schlessinger, J.; Moses, E. Proceedings of the National Academy of Sciences of 
the United States of America 2004 , 101, 3346-3351. 

(13) Paivanranta, B.; Eangner, A.; Kirk, E.; David, C.; Ekinci, Y. Nanotechnology 2011 , 22, 
375302. 

(14) Eirnkes, M.; Pedone, D.; Knezevic, J.; DoILblinger, M.; Rant, U. Nano letters 2010 , 10, 
2162-2167. 

(15) Cooper, D. R.; DaAZAnjou, B.; Ghattamaneni, N.; Harack, B.; Hilke, M.; Horth, A.; Ma¬ 
jlis, N.; Massicotte, M.; Vandsburger, E.; Whiteway, E. International Scholarly Research 
Notices 2012 , 2012. 

(16) Nishiyama, H.; Suga, M.; Ogura, T.; Maruyama, Y; Koizumi, M.; Mio, K.; Kitamura, S.; 
Sato, C. Journal of structural biology 2010 , 169, 438-449. 

(17) Krueger, M.; Berg, S.; Stone, D.; Strelcov, E.; Dikin, D. A.; Kim, J.; Cote, E. J.; Huang, J.; 
Kolmakov, A. ACS nano 2011 , 5, 10047-10054. 

(18) Ei, X.; Cai, W.; An, J.; Kim, S.; Nah, J.; Yang, D.; Piner, R.; Velamakanni, A.; Jung, L; 
Tutuc, E. Science 2009 , 324, 1312-1314. 


14 



(19) Massicotte, M.; Yu, V.; Whiteway, E.; Vatnik, D.; Hilke, M. Nanotechnology 2013 , 24, 
325601. 

(20) Chen, Q.; Smith, J. M.; Park, J.; Kim, K.; Ho, D.; Rasool, H. L; Zettl, A.; Alivisatos, A. P. 
Nano letters 2013 , 13, 4556-4561. 

(21) Yu, V.; Whiteway, E.; Maassen, J.; Hilke, M. Physical Review B 2011 , 84, 205407. 

(22) Jain, P. K.; El-Sayed, I. H.; El-Sayed, M. A. nano today 2007, 2, 18-29. 

(23) Winter, B.; Eaubel, M. Chemical reviews 2006 , 106, 1176-1211. 

(24) Huang, P. Y; Ruiz-Vargas, C. S.; van der Zande, A. M.; Whitney, W. S.; Eevendorf, M. P; 
Kevek, J. W; Garg, S.; Alden, J. S.; Hustedt, C. J.; Zhu, Y. Nature 2011 , 469, 389-392. 

(25) Moghimi, S.; Szebeni, J. Progress in lipid research 2003 , 42, 463-478. 

(26) Upadhyayula, S.; Quinata, T.; Bishop, S.; Gupta, S.; Johnson, N. R.; Bahmani, B.; 
Bozhilov, K.; Stubbs, J.; Jreij, P; Nallagatla, P. Langmuir 2012 , 28, 5059-5069. 

(27) Ein, Y.-C.; Eu, C.-C.; Yeh, C.-H.; Jin, C.; Suenaga, K.; Chiu, P.-W. Nano letters 2011, 12, 
414-419. 

(28) Kochat, V.; Pal, A. N.; Sneha, E.; Sampathkumar, A.; Gairola, A.; Shivashankar, S.; Ragha- 
van, S.; Ghosh, A. Journal of Applied Physics 2011 , 110, 014315. 

(29) De Jonge, N.; Peckys, D.; Kremers, G.; Piston, D. Proceedings of the National Academy of 
Sciences 2009 , 106, 2159-2164. 

(30) Welvaert, M.; Rosseel, Y. PloS one 2013 , 8, e77089. 


15 



